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Abstract 

The escape probability £ x from a site x of a one-dimensional disordered lat- 
tice with trapping is treated as a discrete dynamical evolution by random 
iterations over nonlinear maps parametrized by the right and left jump prob- 
abilities. The invariant measure of the dynamics is found to be a multifractal. 
However the measure becomes uniform over the support when the disorder 
becomes weak for any non-zero trapping probability. 

Implications of our findings in terms of diffusion are discussed. 
05.40+j;05.60.+w;61.43.Hv 
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Diffusion in the presence of disorder and trapping has become by now a classic field. 
This is mainly because the problem per se is mathematically interesting and is well posed; 
also it has a tremendous potential for a wide range of applications which include migrations 
of optical excitations |L|, polymer physics and diffusion-limited binary reactions 0. See 
e.g. [Q for an exhaustive review. 

The standard approach to this class of problems is to write down a second order master 
equation for the probability of the particle to be at a lattice site at a given time, and solve 
it employing analytical or numerical techniques, see for example || . An alternate approach, 
based on the first passage time (FPT) formulation, has attracted growing attention in the 
recent times fl6l]7|]8|,p|,|l0|| . This approach has an advantage in that the master equation is 
first-order to start with. All the transport properties of the system can also be calculated 
from the first passage time formulation. 



Employing FPT formulation for the Sinai model |TT| it was recently shown that the 
distribution of the mean FPT over the disorder exhibits interesting multifractal scaling [|8] . 
Also the probability to escape from one site of the lattice to the next was found to have self 
similar fluctuations P,[IT)|]. The Sinai model however, is an highly idealized, albeit interesting, 
mathematical model, and whose link to physical reality appears to be rather abstract. 

In this work we shall consider a more realistic model for diffusion where a particle diffuses 
by overcoming random barriers but can also be trapped at various sites with site depen- 



dent random probabilities ||12 |. The main characteristic of this class of models is that the 



total probability (called the survival probability) is not a conserved quantity. It has been 
shown that the survival probability is an highly fluctuating function of time and these fluc- 
tuations lead to interesting and unexpected behaviour like enhanced diffusion, breaking of 
self-averaging and emergence of Lifshitz tails [13[ . 



We shall show that this model, when disorder is strong leads to self-similar fluctuations 
of the escape probability, and these can be characterised employing multifractal formalisms. 
However this feature of multifractality disappears when there is no trapping, what ever may 
be the strength of disorder. More importantly, when the strength of disorder goes to zero 



the multifractality disappears even with arbitrary non-zero trapping probability. Purely 
from methodological point of view, we connect diffusion in a trapping environment to an 
Iterated Function System (IFS) ||14|| . Such a formulation, connecting random walks and 
iterated function systems, was proposed very recently in the context of a binary model for 
Sinai disorder |5],|IU|, where we have two maps for random iterations. Here we extend the 
formulation to problems of diffusion on a disordered lattice in the presence of trapping, where 
we have infinity of maps parametrized by the jump probabilities which are chosen randomly 
from a well specified distribution that models the disorder. We restrict our attention to a one- 



dimensional lattice since, as often the case [151, it is contains all the essential characteristics 



of the higher dimensional systems. Furthermore one- dimensional systems are amenable to 
relatively easy analytical and numerical work. 

Let us consider the master equation for the probability G XtX+ i(n) that a particle makes a 
first passage from a site x to a site x + 1 in n steps on a one-dimensional lattice of length N. 
At each site x > 1 we shall indicate by q x 6 [0, 1/2] the probability for making a left jump and 
by p x £ [0, 1/2] the probability of making a right jump (see Fig |IJ). The sojourn probability 
at site x is given by 7(1 — q x — p x ) and the trapping probability is (1 — — q x — p x ). Here 
7 is a parameter which can be continuously tuned from (trapping) to 1 (no trapping) 

The master equation for G x ^ x+ i(n) (x > 1) then reads: 

G XjX+ i{n) = p x 6 1>n + q x G x _ liX+1 (n - 1) + 7(1 - q x - p x )G XtX+1 (n - 1) (1) 

with the boundary condition that site x = — 1 is perfectly reflecting: 

G ,i{ri) = po^i.n + 7(1 - Po)G ,i{n - 1) (2) 

and that G X)X+ i(n) = for x < —1. We assume that {q x , x — 1, N — 1; p x , x = 0, N — 1} 
constitute a set of independent random variables identically distributed in the range (0—1/2), 
and the common distribution is given by 

n(w) = 2 1 " /3 (1 - f3)w- p e(w)9(l/2 - w) (3) 
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where w = q, p and ft G [0, /; 1). Here #(•) is the usual Heaviside function. This distribution 
is known to produce anomalous diffusion when the disorder is strong (ft — > 1) and there is 
no trapping ||. For ft = 0, we find that the distribution is uniform in the range zero to half. 
Thus ft can be tuned from (weak disorder) to 1 (strong disorder). 

Equations (|TJ) and (0) are readly solved employing generating function technique. We 
define 



G^*) = $> n G a ,*+i(n) (4) 

n=0 

Upon the use of convolution theorem, 

Gx-i,x+i{z) — Gx-i,x{z) G X;X +i(z) (5) 

we get: 

G x<x+ i(z) = - — — — - (6) 

1 - 72(1 -q x - p x ) - zq x G x -xA z ) 

for x > 1 and 

This solution has earlier been obtained in Ref. 0. 

We are interested in the behaviour of the escape probability, namely the total probability 
for the first passage from x to x + 1. This is given by: 

+00 

£x = G XjX+ x(z = 1) = 53 G XjX+ x(n) (8) 

n=0 

Using eqns. @ and (0) it is immediately seen that the escape probability satisfies the 
following one-dimensional recursion: 

Cx = Z 7Z X>1 (9) 

1 - 7(1 - q x - p x ) - q x U-i 

with the initial condition 



Eq. (H) can be interpreted as a dynamical map for a fixed p and q. In fact, since p and q are 
random, the evolution £ — ► £1 • • • ~ > £x — > 6r+i —>•••• proceeds by random iteration over 
the maps parametrized by p and q, which are chosen independently and randomly from the 
disorder distribution give by eq. (0) at each stage of iteration. This constitutes an iterated 
function system(IFS), see Barnsley fd| . 

It is immediately seen that when 7 = 1, which corresponds to a lattice with no trapping 
eqns. (§) and fliTf ) lead to £ x = 1, for all x regardless of the choice of {p x , q x }. 

Let us now consider the non-conserved case, for which 7 is less than 1. For given values 
of q and p the fixed point of the map (ffi) is 



[l-7(l-g-p)]~ y/[l-7(l-g-p)] 2 -4^ 
^ 2q ^ J 

and it is stable. It lies (see Fig |) in the region delimited by £* = corresponding to p = 0, 
g > (only left jump) and £* = — 7(1 — p)] corresponding to q = 0, p > (only right 
jumps). 

We now show that the escape probability exhibits self similar fluctuations and these can 



be characterized employing multifractal formalisms [ IB| . From numerical point of view, it 
proves convenient to rescale £ x in such a way that the domain is the interval [0, 1]. To this 
end we employ the standard rescaling: 



£ £min 



£ (12) 



^Max — £min 

We denote by Pi(e), the fraction of the total number of £ values that belong to the i th interval 
of size e = 1/JV. Then the partition function is given by, 

N 

Z(Q,z) = EP?(z) (13) 

i=i 

where the sum is taken over non-empty intervals only. We make the following scaling ansatz, 

Z(Q,e) e ~V^ (14) 
and obtain the scaling exponents t(Q) as: 



t(Q) = lim 



ln[Z(Q,e)] 



(15) 



lne 



Fig. |3] shows a log-log plot of Z(Q, e) versus e(= 1/N) for N ranging from 10 to 3 x 10 6 . The 
linearity of the curves establishes unambiguously the scaling ansatz ([13]). From the scaling 
exponents we calculate the generalized Renyi dimensions, given by D(Q) = t(Q)/(Q — 1). 
Legendre transform of t(Q), defined as 



yields the spectrum of singularities denoted by f(ct). 

Fig. | depicts the scaling exponents t(Q). It is well defined and exhibits clear change in 
slope, establishing that the underlying measure is multifractal. Fig. |5] depicts the spectrum 
of Renyi dimensions D(Q) for various strengths of disorder (/3) and trapping^). We observe 
that when the disorder is strong (/? — > 1), D(Q) remains the same for all values of 7 7^ 1. In 
other words, the strength of trapping does not influence the fractal measures of the escape 
probability, when the disorder in the lattice is strong. However when the disorder is weak 
{(3 — > 0), the spectrum of Renyi dimensions changes from one trapping rate to the other. 
Also in the limit of /3 — > 0, the D(Q) curve becomes flat with unit intercept, for all values 
of 7 7^ 1, implying that the measure is uniform and space filling. To capture in a simple 
fashion the dependence of the fractal measure on the strength of disorder, we depict in Fig. 
^ the variation of the information dimension D(l) and the correlation dimension D(2), as 
a function of (3 for a fixed value of 7 = 0. We find that both D(l) and D(2) decrease with 
increasing strengths of disorder. We plot in Fig. [7| the / — a curve. It is worthwhile noticing 
that since the slope of t(Q) for Q — > —00 saturates at unity, the side of f(a) for a > 1 does 
not exhist. 

A natural question that arises in this context relates to the implications of our findings 
to the transport properties of the disordered systems. More specifically we ask the question: 
Is there a connection between anomalous diffusion and fractal fluctuations? For example 



f(a) = aQ- t(Q) 



(16a) 



(16b) 
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it has already been shown analytically that even in the absence of trapping, the system 
exhibits anomalous diffusion when the disorder is strong (/3 — » 1). However when the 
disorder is weak, the anomaly in the diffusion process disappears. On the other hand, our 
finding is that when the disorder is strong, for arbitrary trapping rate, (so long it is non 
zero), the escape probability exhibits multifractal fluctuations. However, when the disorder 
becomes weak {(3 — > 0), the multifractal features disappear. More studies with different 
disorder models are required to shed more light on the issues raised here. We hope this work 
would spur some activities along these directions. In any case we believe, it would worth 
the effort to investigate the intriguing connection between anomalous diffusion and fractal 
fluctuations in random and trapping environments. 
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FIGURES 

FIG. 1. Definition of the hopping probabilities. At site x the left jump probability is q x and 
the right jump probability is p x . The sojourn probability is S x = 7(1 — q x — p x ) while the trapping 
probability is T x = (1 - 7)(1 - q x - p x ) where 7 € [0, 1]. 

FIG. 2. The map £ x = for P = 0.2, q = 0.1, p = q = 0.15 and p = 0.1, g = 0.2. Fixed 

points are the interections of the maps with the £ x = £ x -\ line. 

FIG. 3. Behaviour of the partition function Z(Q, e) vs. e in a log-log plot for (3 = 0.3 and 
7 = 0.99. Q varies are from —5 (top curve) to +10 (bottom curve) in units of 1 with Q / 1. 

FIG. 4. Plot of t(Q) vs. Q. 

FIG. 5. Spectrum of the Renyi dimensions D(Q) vs. Q. 

FIG. 6. The information dimension D(l) and of the correlation dimension D(2) as a function 
of the strength of the disorder (3, for 7 = 0. 

FIG. 7. The spectrum of singularities for (3 = 0.7, 7 arbitrary and (3 = 0.3, 7 = 0. 
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